Revisiting the space-time gradient method: A time-clocking perspective,high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient (STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric Gauss-Seidel/Gauss-Seidel (LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-count-reduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.     

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

GNSS驯服芯片级原子钟方法研究CSCD

芯片级原子钟是一种体积小且功耗低的高精度时钟源,具有广泛的用途。针对这一特点,设计了基于GNSS的芯片级原子钟驾驭算法。以GNSS系统时作为参考,测量芯片级原子钟与GNSS系统时间的钟差,并对芯片级原子钟进行钟差建模,获取其特征参数。通过乒乓法计算出钟驾驭调整量,对芯片级原子钟进行控制,最终将芯片级原子钟驾驭到GNSS系统时间上。经过实验验证,在驾驭时间常数为100s的情况下,芯片级原子钟与GNSS系统时间的时钟同步误差在-7.5~7.5ns之间;1h频率准确度为5.8×10-13;平均时间为10000s时的频率稳定度为3×10-13。     

基于重叠网格的唇口移动变几何进气道数值模拟研究

吸气式冲压发动机的进气道需要在比较宽广的马赫数范围内工作,而变几何进气道具有良好的气动性能,成为了最佳选择。重叠网格技术是解决存在复杂运动边界流动模拟问题的有效策略,在应用于变几何进气道非定常流动模拟时需要进一步完善。采用发展的重叠网格技术,对变几何平动方式的进气道开展了数值模拟研究。针对存在物面接触重叠网格的挖洞操作在执行时无法取得合理的插值模板问题,引入了虚拟单元赋值方法和“Collar”网格方法,实现了多体相交时重叠网格的计算。为了引入“Collar”型重叠网格,对其算法作了适当修改以适应本论文所建立的洞面优化方法。对平移过程中内收缩变化所造成的对启动过程的影响进行了研究。通过迎着来流方向平移唇口,增大了进气道的内收缩比,导致内压缩段出现边界层分离,分离泡随着唇口的前移而增大,当其移出内压缩段时,原本启动的进气道进入不启动状态,而后移唇口减小了内收缩比,使分离泡减小并向下游移动,辅以边界层排移,可以改善进气道的启动性能。此外,平移速度和振荡频率对启动过程的影响也进行了研究,进气道迟滞效应随平移速度和振荡频率增大而增大。     

后摆心潜入式喷管柔性接头力学特性数值仿真

针对后摆心潜入式喷管柔性接头在燃气压力与摆动联合载荷下,弹性件与增强件的力学特性计算问题,建立了对称刚体-柔性体结合的几何模型;采用了有限元方法,获取了柔性接头在燃气压力、摆动载荷及联合载荷下应力、弹性比力矩分布规律。研究表明:随着燃气压力增大,弹性比力矩降低;随摆角增大,弹性比力矩增大。文中的数值仿真方法,可以应用于柔性接头力学特性的计算,为柔性接头强度的安全性评估提供指导。     

作动器故障下的无人机容错控制方法

在执行任务过程中,无人机的传感器、作动器等均可能出现故障。文章针对常规布局无人机的作动器故障,提出了 1种反步法和控制分配相结合的容错控制方法。首先,建立无人机数学模型,并对作动器故障进行分类和建模;然后,根据模型设计反步最优控制器和基于控制分配的容错控制器;最后,通过仿真验证表明,所设计的容错控制方法能够实现作动器故障下的姿态快速稳定控制,且稳定性好,基本无侧滑角,各操纵面均在约束范围内,达到容错控制要求。     

自由转子陀螺仪的伺服测试技术研究

自由转子陀螺仪是目前精度最高的一种陀螺仪,没有精密的力矩器,因此,只能采用双轴伺服法来辨识其漂移误差模型。介绍了双轴伺服测试法的基本原理,给出了双轴伺服转台的运动轨迹及重力矢量的变化规律。根据转台的运动特性,建立了自由转子陀螺仪的漂移误差模型,导出双轴伺服转台转角数据估计漂移误差模型系数的算法。从测试结果可以看出自由转子陀螺仪的长时间精度高,表明采用伺服测试法能够获得极高的测试精度,伺服转台可以作为测试自由转子陀螺仪精度的装置。     

基于双目视觉的返回舱着陆横向速度测量

通过分析返回舱在着陆时所具有的水平速度也可能对着陆安全构成威胁的事实，论述了为返回舱增加横向缓冲装置的必要性，由此提出亟需解决的问题——在返回舱着陆下降过程中测量其水平运动速度。提出了全新概念的以视觉图像为唯一信息来源的返回舱水平速度测量方法，仅以双目摄像机系统作为传感器，避免了对返回舱原有传感器系统的改造。通过对着陆区域地面图像的采集和在线处理，利用双目交会测量原理得到返回舱当前的大地坐标和姿态。仿真实验证明，方法简单易行，测量结果可靠。     

固体发动机药柱粘结试件的三维应力分析

推进剂／衬层／绝热层矩形粘结试件已作为随发动机测试试件，用于监控药柱最薄弱的推进剂／衬层界面粘结质量。由于试件较厚，又材料具有粘弹特性，因此应对它进行三维粘弹性有限元分析。本文分析了推进剂／衬层界面附近的应力分布情况，并给出试件启裂点有效应力与拉伸平均应力之比的集中系数，供药柱结构完整性分析人员使用     

定速拉伸的固体推进剂双轴板条的应力分析

对固体推进剂双轴板条试件进行了有限元分析，并考虑了粘弹性对应力结果的影响，分析结果表明，双轴板条试件中心域为平面应变状态，与发动机药柱中部内孔表面相似，并给出了计算杨氏模量所需的修正系数。